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In the framework of continuum theory we study orientational transitions induced by electric and magnetic fields in ferronematics, 
i.e., in liquid-crystalline suspensions of ferromagnetic particles. We have shown that in a certain electric field range the magnetic 
field can induce a sequence of re-entrant orientational transitions in ferronematic layer: nonuniform phase — uniform phase — 
nonuniform phase. This phenomenon is caused by the interplay between the dipole (ferromagnetic) and quadrupole (dielectric 
and diamagnetic) mechanisms of the field influence on a ferronematic structure. We have found that these re-entrant Freed ericksz 
transitions exhibit tricritical behavior, i.e., they can be of the first or the second order The character of the transitions depends 
on a degree of redistribution of magnetic admixture in the sample exposed to uniform magnetic field (magnetic segregation). We 
demonstrate how electric and magnetic fields can change the order of orientational transitions in ferronematics. We show that 
electric Freedericksz transitions in ferronematics subjected to magnetic field have no re-entrant nature. Tricritical segregation 
parameters for the transitions induced by electric or magnetic fields are obtained analytically. We demonstrate the re-entrant 
behavior of ferronematic by numerical simulations of the magnetization and optical phase lag. 
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1 Introduction 

■It is known, that liquid ferromagnetics do not occur in nature 
'since melting temperature of all known substances is always 
higher than the so-called Curie temperature at which the long- 
range order in orientation of magnetic moments of atoms is 
destroyed and the phase transition to paramagnetic state takes 
place. For this reason liquid ferromagnetics can only be the 
artificially fabricated media. First were synthesized the so- 
called magnetic fluids ' which are colloidal suspensions of 
magnetic particles in isotropic liquids. These suspensions do 
not possess ferromagnetic properties since after the external 
magnetic field is off they do not retain the residual magnetiza- 
tion. Magnetic fluids are inherently paramagnetic and because 
of their high magnetic susceptibility are often referred to as 
"superparamagnetics". 

However, it is possible to synthesize a suspension with fer- 
romagnetic properties on the basis of an anisotropic liquid. 
'One of such materials is ferronematic^'^ (FN), i.e., low con- 
centrated suspension of anisometric magnetic particles in ne- 
matic liquid crystal (NLC). In order to create a noncompen- 
sated ferronematic possessing spontaneous magnetization, we 
need to add anisometric magnetic particles into the liquid crys- 
tal heated above the clearing point and then cool it to liquid 
crystalline state in an external magnetic field. As a result, the 
ferronematic will possess magnetization even in the absence 
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of a magnetic field. Due to the liquid crystal matrix FN has 
good fluidity and anisotropy of physical properties, and em- 
bedded in NLC-matrix ferroparticles stipulate a strong mag- 
netic response of the suspension. 

In addition to fluidity and high sensitivity to the external 
fields, ferronematics demonstrate the phenomenon of mag- 
netic admixture redistribution in a uniform magnetic field. 
Since the magnetic particles are not fixed to the NLC-matrix, 
they have the opportunity of spatial movement migrating to 
areas where the sum of their magnetic and orientation ener- 
gies in liquid crystal matrix is minimal. The phenomenon of 
magnetic particles accumulation in "favorable" parts of the 
sample under a uniform magnetic field is called the segre- 
gation effect.^ This effect leads to tricritical behavior both 
ferrocholesteric-ferronematic transition and magnetic Freed- 
ericksz transition in ferronematics.^ 

Due to coupling of magnetic particles with liquid crystal 
molecules there exist two mechanisms of magnetic field ef- 
fect on the orientational structure of the suspension: the mag- 
netic quadrupole (influence on the liquid-crystal matrix) and 
magnetic dipole (influence on the magnetic particles). The 
competition between these mechanisms leads to the thresh- 
old Freedericksz-like transition in unbounded ferronematics-i 
and re-entrant orientational transitions in ferrocholesterics.— 
An electric field directly acting only on the LC-matrix gives 
one more mechanism of influence — the electric quadrupole. 
The competition between these three mechanisms additionally 
leads to re-entrant orientational transitions in ferronematics.— 
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As is shown in Ref. the re-entrant transitions are possible 
in ferronematics subjected to the magnetic field without elec- 
tric field, but under special bistable coupling conditions on the 
bounding plates. 

In recent years there have been many experimental ' '"'^ and 
theoretical ^^"'^ works on the physical properties of magnetic 
suspensions based on liquid crystals and the peculiarities of 
their behavior in external fields. A detailed review of the 
history of research, classification, the most interesting exper- 
imental and theoretical results, as well as numerous applica- 
tions of liquid-crystalline nanocolloids can be found in Ref.—. 

In this paper we study magnetic- and electric-field-induced 
re-entrant orientational transitions in ferronematics. Taking 
into account segregation effects leading to the tricritical be- 
havior,^ we emphasize on interplay between dipolar (ferro- 
magnetic) and quadrupole (dielectric and diamagnetic) modes 
of field influence in ferronematics. Our aim is to study the role 
of these interactions as well as segregation effects in tricritical 
behavior of re-entrant orientational transitions in ferronemat- 
ics. We ask ourselves: is it possible to change the order of the 
Freedericksz transition in ferronematics by external fields? 

The rest of the paper is organized as follows. Section [^in- 
troduces a set of equations describing the orientational struc- 
ture of ferronematic in electric and magnetic fields taking into 
account the segregation effect. In Sec. |3]we analyze the re- 
entrant magnetic Freedericksz transitions in electric field and 
electric transitions in magnetic field using Landau expansion 
of ferronematic free energy. Section 0] is devoted to orienta- 
tional, magnetic and optical properties of ferronematics. Con- 
cluding remarks are presented in Sec. |5] 



2 Orientational equilibrium equations 

We consider a layer of ferronematic liquid crystal having 
thickness of L, placed between two parallel plates (see Fig. 
Let us introduce a rectangular coordinate system, axis x is di- 
rected along the plates, axis z — perpendicular to the plates; 
we choose the origin in the middle of the layer CoupUng 
of the director n to the layer limiting plates we assume to be 
rigid (i.e., orientation of the director at the boundary does not 
change under the influence of the external field) and planar 
(axis of easy orientation on the layer surface no is parallel 
to the axis x). We apply electric £ — (0, 0, £) and magnetic 
= (0, 0, M) fields orthogonally to the ferronematic layer. 

Equilibrium state corresponds to the minimum of free en- 
ergy functional 3^ ~ J FydV. The bulk free-energy density 
Fy of ferronematic in electric and magnetic fields under soft 
surface coupling of magnetic particles to NLC-matrix can be 




Fig. 1 Orientation of ferronematic layer in electric £ and magnetic 
M fields 



written as foUows-^i^ 
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Fi 
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F6 



F1+F2 + F3+F4 + F5+F6, 
1 



(1) 



[Ki (V ■nf + K2{n-Vx n)^ + Ki{nxVx nf] 



-ieoe„(n-£)^ 



F3 = -^oM,/m-JC, 



F,= '-f{n-m)\ 
a 



Here Ki , K2 and Kj, are elastic modules of NLC (Frank con- 
stants), Xa is the anisotropy of diamagnetic susceptibility of 
NLC, /io is the permeability of vacuum, is the saturation 
magnetization of magnetic particles material, eo is the permit- 
tivity of vacuum, e„ is the anisotropy of LC permittivity, /(r) 
is the local volume fraction of magnetic particles in the sus- 
pension, m is the unit vector of suspension magnetization, v 
is the volume of a ferroparticle, d is the diameter of a ferropar- 
ticle, kB is the Boltzmann constant, T is the temperature, w is 
the surface density of coupling energy between nematic LC 
molecules and the surface of magnetic particles. 

The term Fi in expression ([U represents bulk free-energy 
density of orientational elastic deformations of the director 
field (Oseen-Frank potential), F2 is the bulk free-energy den- 
sity of the magnetic field interaction with nematic matrix (it 
describes the quadrupole mechanism of the magnetic field in- 
fluence on FN), is the bulk free-energy density of the mag- 
netic field interaction with magnetic moments of ferroparticles 
(dipole mechanism of the magnetic field effect on FN), F4 is 
the contribution of the mixing entropy of the ideal solution of 
magnetic particles to the bulk free-energy density, F5 is the 
bulk density of the surface coupling energy of magnetic parti- 
cles with the director, is the bulk free-energy density of the 
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interaction of electric field with the nematic. 

We assume a surface density of coupling energy w > 0, so 
that in the absence of external fields the minimum of free en- 
ergy corresponds to mutual orthogonal orientation of the di- 
rector and magnetization (n_Lm). Let us consider ferrone- 
matic with positive anisotropics of diamagnetic susceptibility 
Xa and dielectric permittivity £„. Because of initial orthogonal 
orientation of the director and magnetization there appears a 
competition between dipole magnetic fioMsfm ■ 5t] and 
quadrapole magnetic l^oXai^ ■ ^)^] as well as electric 
eoea(n • £)^] mechanisms of influence in the presence of 
electric and magnetic fields. Assuming a small volume frac- 
tion of ferroparticles in ferronematic (/ ^ 10^^), but still suf- 
ficient for collective behavior of the suspension 22,23^ magnetic 
dipole-dipole interactions can be neglected. 

In the geometry under consideration (Fig. [1} the director n 
and unit vector of magnetization m can be written as 



n = [cos(p(z), 0, sin(p(z)], 
[-sinv/(z),0, cosv/(z)]. 



(2) 



To analyze the problem it is convenient to work with dimen- 
sionless variables. We choose layer thickness L as a unit of 
length and introduce the dimensionless coordinate z = z/L. 
Further on, for simplicity, we omit the tilde sign over the di- 
mensionless variables in the equations. Let us define the di- 
mensionless fields 



H = 'KL 



£()£a 



(3) 



and dimensionless parameters 



b = MJL 



Mo 

XaKl 



Kid 



kBTjL^ 



(4) 



Here / — Nv /V is the average volume fraction of magnetic 
particles in FN. As a unit of magnetic field 5{ in (l3]l we choose 
value = L~^^Ki/{iM)Xa)- In this case at « the en- 
ergy of elastic distortions Fi (Oseen-Frank potential) and dia- 
magnetic term F2 in the free energy of ferronematic dl) appear 
to be of the same order. For J-C > orientational distortions 
are due to diamagnetic anisotropy of LC-matrix (quadrupole 
mechanism of magnetic field effect on FN). A similar com- 
parison of the elastic Fi and dipole ^3 contributions gives an- 
other characteristic value of the field = K\/{pL()MsfL}). 
For 'K > orientational distortions are due to the influence 
of the magnetic field on the particles (dipole mechanism). Pa- 
rameter b — 'Kq/'Kj is the ratio of these characteristic fields 
and characterizes the mode of magnetic field effect on FN.— 
For b> I {'Kq > 3-C11) distortion of FN orientational structure 



in weak fields is caused by the dipole mechanism, whereas 
for b < I {'Kcj < Jid) — quadrupole mechanism. Change of 
the mode of influence from dipole to quadrupole (and vice 
versa) takes place in the fields for which the contributions 
F2 and Fj, into the free energy are of the same order, i.e., at 

0<^:Ko^Mj/Xa- 

As a unit of electric field £ in (|3]l we have chosen the value 
£^ = L^' y^Ki/{£()£a), which gives the characteristic value of 
Freedericksz transition field in nematic. It is determined from 
the comparison of elastic Fi and dielectric contributions 
into free energy ([T). 

We have defined the so-called segregation parameter— 
K ~ (L/A)^ which is the square of the ratio of two char- 
acteristic lengths: layer thickness L and segregation length 

A = {vKi/ksTf) '-. The last one characterizes the magni- 
tude of concentrational stratification area in FN. For >ir ^ 1 
the distribution of magnetic particles in a ferronematic layer 
is close to uniform because characteristic dimensions of the 
area where concentrational redistribution takes place becomes 
small in comparison with the thickness of the layer For < 1 
the nonuniformity of magnetic particles distribution in a layer 
becomes significant. Moreover, we introduce dimensionless 
coupling energy of magnetic particles with the NLC-matrix a 
and the anisotropy coefficient of orientational elasticity k. 

It has been already noted, that FN has two mechanisms of 
response to the applied magnetic field. The first one (F3), lin- 
ear to the field H, stipulates the orientational behavior of FN 
in weak magnetic fields. This contribution describes the influ- 
ence of an external magnetic field directly on magnetic mo- 
ments of ferroparticles, and due to coupling F5 indirectly on 
the nematic matrix. The second mechanism {F2), quadratic to 
the field H, corresponds to the influence on a diamagnetic LC- 
matrix and, via it, in accordance with F5, — to the magnetic 
moments of ferroparticles. However, in some cases it is neces- 
sary to separate the influence onto magnetic particles and onto 
the LC-matrix. This can be done by placing FN into the elec- 
tric field, the influence of which on LC-matrix is described by 
term fg. 

Orientational part of the full free energy of a ferronematic 
3^ ([T]l in dimensionless form can be written as 



F = -L:,. 

KiS 

1 



1/2 
-1/2 
2 , 



bHg cos i/A- 



-{H^+E^) sin^ (p + >cglng + ag sir? {(p-w) 



dz, (5) 



where S is the surface area of the plates limiting FN layer, also 
we introduced functions 



g{z)^f{z)/f, 
3C(^) ~ cos^^ + A:sin^^. 



(6) 
(7) 
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Here g{z) is the reduced volume fraction of magnetic particles 
in FN. 

Free energy (|5]l is a functional over functions (p{z), y/{z) 
and g{z)- Minimization of the free energy (|5} with respect to 
(p{z) and \j/{z) gives equations determining orientation angles 
of the director and magnetization: 



-E^) sin2(p 



-ff§sin2(^- V/) = 0, 
fe//sin V/- (7sin2(9 - i/) = 0, 



with the boundary conditions 



(p(-l/2) = (p(l/2)=0, 



(8) 
(9) 



(10) 



corresponding to the rigid planar coupling of the director on 



the surfaces n 



=±1/2 



(1,0,0). Hereafter the prime denotes 



the derivative with respect to the dimensionless coordinate z- 
Equilibrium distribution of magnetic particles over the FN 
layer is described by minimising the free energy functional (|5]) 
with respect to g{z) under the condition of constant number of 
particles in a suspension / fdV = Nv: 



(bH 



g = gexp<{ — cosv/ sin^(^-v/)^, (11) 



where 



Q^= / exp< — cosy/ sin^ {(p - > dz. (12) 

J-l/2 L ^ ^ J 

Note that the so-called couphng equation (|9|l determines 
mutual orientation of the director and magnetization,— 
whereas eqn ( fTTT i describes the segregation effect,^ consist- 
ing in magnetic particles concentration growth in those parts 
of the sample where the sum of their magnetic energy in the 
field K and orientational energy in LC-matrix is minimal. For 
1 magnetic segregation effects can be neglected. In this 
case it is obvious from eqns ( fTTT i and (fT2] l that normalization 
integral I, whereas volume fraction of magnetic particles 
g(z)->l,/.e.,/(z)^7. 

Let us make estimates of the dimensionless variables (HJl 
using typical values of material parameters for nematic liq- 
uid crystals ^'^^ and magnetic particles^: Xa =2.1 x 10^^, 
Ki = 6.4 X lO-'^N, K3 = 1.0 X lO-'^N, T = 298K, / = 
2.0 X 10-\ M, 5 X lO^Am-i, w = lO"*^ - lO^'^Nm-', 
d = 1.5x lO^^m, V = 8.8 X lO^^^m^, and taking layer thick- 
ness L=250^m, we find/tsa 1, 1, CJ « 10"^ - 1, « 10 
and K fa 10^^. As can be seen from these estimates, the small- 
ness of K testifies to the importance of segregation effects in 
the problem under consideration. 

Equations ([8]l, Q, and ( fTTT ) with boundary conditions ( fTOl i 
admits a uniform solution (p{z) = yf{z) = corresponding to 



the planar (n = no) texture of FN with orthogonal (n_Lin) ori- 
entation of the director and magnetization. At the same time 
they determine also nonuniform orientations of the director 
and magnetization of FN under mutual influence of electric 
and magnetic fields. Let us consider nonuniform solutions for 
the orientational angles of the director and magnetization. At 
first we multiply eqn (|8]l by (p' and eqn (|9|l by g\j/', then we 
subtract the second equation from the first. As a result we 
obtain 



d 
dz 



Q. (13) 



Taking into consideration that in the middle of the layer the 
deviation of the director from the axis of easy orientation is 
maximal, i.e., (jo' = at z = 0, the first integral of eqn (fT?t 
takes the form 

X{(p){(pf^ {H^ +E^) (cos^ (p - cos2 (po)+2>c {go - g) . 

(14) 

Here go = g{(po,nfo) is the distribution function of magnetic 
particles, ^0 = ^(0) and \f/o = (/(O) are angles of director and 
magnetization orientation in the middle of the layer, respec- 
tively. 

Integration of eqn (fT4l i for z > with boundary conditions 
( [Tot gives 



1 



where 



X{(p,xff) = 



y/5l{(p,\l/)d(p = --Z, 

2 



cos^ (p + k sin^ (p 



(15) 



(//2 + £2) (cos2 (p - C0S2 (Pq) + 2>C {go - g)' 

(16) 

We choose plus sign in eqn (flST l, that corresponds to the posi- 
tive (counterclockwise rotation) values of the angle of the di- 
rector orientation. 

In the middle of the layer (z = 0) the angle (p has maximum 
value (po, which is determined by the following equation, aris- 
ing from (flST l 

' ^3i{(p,xi/)d(p^^. (17) 



Transforming the normalization integral Q (fT2] l in eqn (fTTT i 
with the help of expression (fT4l i. we get the equation for the 
distribution function g{z), 



f<PO 

/ g{(p,\lf)V^i<P,¥)d(p 
Jq 



(18) 



Thus, eqns ([TSb-dTSll and coupling equation (|9]l determine the 
angles of orientation of the director (p and magnetization xj/, 
and distribution function g{z) = f{z)/f of particles as func- 
tions of magnetic H and electric E fields and material param- 
eters of ferronematic. 
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3 Reentrant Freedericksz transitions in fer- 
ronematics 



The director rotation angle is defined by minimizing of free 
energy (ISTT l with respect to the 



Equations dHJ, (|9]l, and ( fTTT i with boundary conditions (fTOl l ad- 
mit a uniform solution (p{z) = yf{z) — 0, corresponding to the 
planar (n — no) texture of FN with orthogonal (n±m) orien- 
tation of the director and magnetization. However, this solu- 
tion becomes unstable and is replaced by nonuniform solution 
when the external field (electric or magnetic) reaches a thresh- 
old value. Such transition from uniform to nonuniform state is 
called the Freedericksz transition^. It is the result of competi- 
tion of surface forces that orient the liquid crystal on the layer 
boundary, forces of orientational elasticity and external force 
fields. In the considered geometry (Fig. [T) the Freedericksz 
transition in FN layer can be caused by both magnetic and 
electric fields. Let us consider first the magnetic-field induced 
nonuniform orientation of the director at a given electric field. 



3.1 Magnetic Freedericksz transitions in an electric field 

Close to the Freedericksz transition, angles (p{z) and \j/{z) are 
small, so in the lowest order expansion of eqns (jSJ, (|9]l, and 
dTTb with boundary conditions ( fTOl i we find 

(p{z) = (po cos 7ZZ, \ir{z) =SH(p{z), SH = 2a/{2a + bHc), 

(19) 

where 90 *C 1, and critical magnetic field is determined by 
equation^ 



^2 



72 



SHhH,. 



(20) 



At zero magnetic field, // = 0, it gives the threshold value of 
the electric field Ec in FN which coincides with the Freeder- 
icksz field E^'- = ;r in pure nematic. ^ At zero electric field, 
£ = 0, eqn ( |20] | reduces to the expression obtained in Ref. ^. 

Free energy (|5) can be expanded as a power series in the 
(p{z) and ( fT9] l. The fourth-order expansion in the Landau 
form gives 



(21) 



where 



Fq = -bH, «// = — [H^{2-SH)+SH{n^~E^)] , 



Here we introduce 



[H^+E^-n^y 



(22) 



(po = ±JrH 



H H, SK[H^{2-SH)+SHin^~E^)] 

Yh — 

x->^H 3b[HcSH{2-SH)f+^TZ^kH, 

(23) 

As seen from formula (|23), segregation parameter corre- 
sponds to the tricritical point, in which the character of tran- 
sition between the uniform and nonuniform phases changes 
from the second order to the first. For k > (weak segrega- 
tion) at H — He orientational transition is the second order, 
whereas for k < x^j (strong segregation) the transition be- 
comes the first order In the latter case, ferronematic exhibits 
orientational bistability. At zero electric field, £ = 0, formulas 
i2l[ - ( |23] | coincide with expressions, obtained in Ref.^. 



nonuniform 




1 2 3 4 5 

magnetic field, H 

Fig. 2 Bifurcation phase diagram of orientational transitions in FN 
for a = 10 and b = 10. Here Ej^'- = Hj/' = n am fields of electric 
and magnetic Freedericksz transitions in pure nematic, 77 is the 
width of the region where the re-entrant phases are possible 

The bifurcation phase diagram (|20] i of orientational transi- 
tions in FN is plotted in Fig.|2] In the case of weak magnetic 
segregation {x > x^j) the curve E{H) defines the boundary 
of the second-order Freedericksz transitions in FN induced by 
the joint action of electric and magnetic fields. The region 
under the curve E{H) corresponds to the unperturbed state of 
FN, i.e. , uniform planar (n = no) texture of FN with mutual or- 
thogonal orientation of the director and magnetization (n±ni), 
whereas the region over the curve corresponds to the nonuni- 
form state. 

Characteristic feature of the phase diagram (Fig. |2) is the 
presence of re-entrant phase transitions in FN (nonuniform 
phase — uniform phase — nonuniform phase) caused by com- 
petition between quadmpole electric £oe„ (n • £)^] and mag- 
netic /ioXa(n IK)^] mechanisms of influence on NLC- 
matrix, and dipole ^M^fva ■ 3<] influence of magnetic 
field on magnetic moments of ferroparticles. As can be seen 
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from Fig. [2] for E^r'^ < E < E^ + Tj at i/ = ferronematic is 
in nonuniform phase because E^ is a field of electric Freed- 
ericksz transition. Application of a magnetic field induces the 
transition to uniform phase and then with further increase ofH 
uniform phase again is changed to nonuniform. Thus, critical 
field He for a given E has two values: Hi at the transition from 
nonuniform state to uniform, and H2 > Hi — at the re-entrant 
transition from uniform state to nonuniform (Fig. |5^). In the 
first case, formula (l2Jt has parameter Jh < 0, in the second 
— pai-ameter > 0. For E > E^ + rj Freedericksz transi- 
tion can not be induced by a magnetic field. The width of the 
region where the re-entrant phases are possible 



ri = 



- . 77-2 



2(7)4 la{q-2af 



6^(7+ (^-2(7)2 



(24) 



where 



8 C7 + 27^2 _^ 3^7 J 3 ( 1 6 C7 + 27/72) 



was obtained analytically in Ref. ^. As is seen from (l24l i the 
interval of electric field values, admitting re-entrant uniform 
phase of FN, expands with the increase of coupling energy 
(7 and the parameter h, and does not depend on segregation 
degree, which is characterized by the parameter >f. 

For the magnetic Freedericksz transition in a non-zero elec- 
tric field, the tricritical value >r^ of segregation parameter (l22t 
is a function of £. As the electric field increases the param- 
eter >c*^ decreases as shown in Fig. |3] Thus at zero electric 
field the magnetic Freedericksz transition in FN is the first- 
order transition {k < x'^). The application of electric field 
can decrease the value so that will be satisfied the inequal- 
ity >f > and magnetic Freedericksz transition becomes of 
the second order Such behavior of the orientational transi- 
tion is confirmed by numerical simulation of eqns (|9), (flST l 
- (fTSl l and shown in Fig. Moreover, on phase diagram 
(Fig. |2]i there are regions of ambiguity, in which fixed elec- 
tric field corresponds to different values of >f^. Decrease of 
orientational elasticity anisotropy k leads to the increase of tri- 
critical value of segregation parameter, besides, the width 
of ambiguity region Tj does not depend on anisotropy con- 
stant k (Fig.[3t) which is consistent with expression (|24] |. De- 
crease of parameter b leads to the decrease of and the 
narrowing of ambiguity region (Fig. Isji). The latter disap- 
pears only at = 0, that corresponds either to pure NLC, or to 
non-magnetic admixture. For parameters marked at the plane 
{k,E) in Fig. [3b by points Mi = (0,1.2), M2 = (4.0,1.2), 
Pi = (3.8,2.0), P2 = (3.8,0.4), P3 = (3.8,0.25), the director 
orientation angle (po in the middle of the layer is shown in 
Fig.lSl 
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Fig. 3 Tricritical value of segregation parameter as a function of 
electric field E for different values of (a) orientational elasticity 
anisotropy constant k and (b) parameter b, which characterizes the 
modes of magnetic field influence on ferronematic; Mi = (0, 1.2), 
M2 = (4.0, 1.2), Pi = (3.8,2.0), P2 = (3.8,0.4), P^ = (3.8,0.25) 



3.2 Electric Freedericksz transitions in a magnetic field 

Let us now consider the Freedericksz transitions between uni- 
form and nonuniform phases induced by electric field in the 
presence of a magnetic field. From expression (|20i and Fig. |2] 
it follows that the threshold electric field Ec of Freedericksz 
transition in FN is a single-valued function of H: 



Ec = y/^^ + bHsE-H^, se ^ lo / {2a + bH), (25) 

so, the re-entrant orientational transitions can not be induced 
by the electric field. Moreover, magnetic field H should not 
exceed the critical field He of the magnetic Freedericksz tran- 
sition ^ in FN at £■ = 0, otherwise the Freedericksz transition in 
FN can not be induced by electric field (Fig.|2]i. If this require- 
ment is satisfied, free energy (|5} close to E^ can be represented 
as a power series in <p{z) — (pocos(7rz) and (/(z) = se<p{z), 
where ^0 ^ 1 • The fourth-order Landau expansion of the free 
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energy takes the form 



■ Fq - 



as 



PE 4 



(26) 



where 



-hH 
1 

16^ 

Here we introduce 



Fq 



{bHsEY 
3bHsl {2~SEf+47!:^k' 



(27) 



Minimization of free energy (1261 1 with respect to gives the 
expression for the director rotation angle in the middle of the 
layer 



(po^±\ Ye 



E-E, 



\6xEc 

7£ = r7 T- > 0- 

3bHsl {2-se) +47t^k 



(28) 

Formula (|28] l shows that electric Freedericksz transition in a 
magnetic field also demonstrates tricritical behavior For x > 
(weak segregation) orientational transition is the second- 
order transition, and for x < >t*^ (strong segregation) — the 
transition of the first order The difference of Freedericksz 
transitions in FN induced by the electric field in the presence 
of magnetic field from magnetic Freedericksz transitions in 
electric field is that these transitions can not be re-entrant tran- 
sitions, because critical field Ec is a single-value function ( |25] ) 
of a magnetic field H. At the zero magnetic field, H = Q, 
electric orientational transitions in FN are the second-order,— 
because >r|- = 0. 

Tricritical value of the segregation parameter (l27t for 
the electric Freedericksz transition in FN (1281 depends on the 
magnetic field H. As the magnetic field increases the tricritical 
segregation parameter increases as shown in Fig. S] Thus, if 
at zero magnetic field the electric Freedericksz transition in 
FN is a second-order transition (i.e., x > k'^), the magnetic 
field will lead to the increase of k'^ so that x < and the 
electric Freedericksz transition in a magnetic field will become 
the first order (Fig. |6). 

Moreover, decrease of orientational elasticity anisotropy 
k increases the tricritical value of segregation parameter >r|- 
(Fig. 2^), whereas decrease of parameter b, per contra, leads 
to the decrease of (Fig.|4}3). 

3.3 The field of the first-order equiUbrium Freedericksz 
transition 

As stated above, magnetic orientational transition between 
uniform and nonuniform states of FN in the electric field for 
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Fig. 4 Tricritical value >c'^ of the segregation parameter as a 
function of magnetic field H for different values of (a) orientational 
elasticity anisotropy constant k and (b) parameter b, which 
characterizes the modes of magnetic field influence on ferronematic 



X < x^ or electric transition in a magnetic field for x < x"^ 
are the transitions of the first order Therefore equilibrium 
transition fields is not the same as E^ or H^. Let us determine 
magnetic //, and electric E, fields of equilibrium Freedericksz 
transitions of the first order These fields are determined from 
the equality of the free energies of perturbed (|5} and unper- 
turbed states of ferronematic 

F=Fq, (29) 

where the free energy of unperturbed state 

Fo = F\^=^=Q = -bH. (30) 

Using eqn ( fT4b expression ( |29] l is reduced to 

Jo 2 (31) 



-x[lneQ-gQ]+bH = 0, 
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where e is the base of the natural logarithm, whereas functions 
3<!((p), I/), i/), and normalization integral Q are de- 
termined by expressions (jTjl, (fTSI l, (fT¥t . and (fTZt respectively. 
This equation is solved together with the equations of orien- 
tational equilibrium (|9]l, ( fTTT i. and (fTST l. For magnetic Freed- 
ericksz transition in a given electric field E eqn dSTT i allows 
to obtain //, (Fig. |5]l, and in the case of electric Freedericksz 
transition at fixed magnetic field H it determines £, (Fig.|6]l. 

4 Orientational, magnetic and optical proper- 
ties of ferronematic 

4.1 Magnetic Freedericksz transitions in an electric field 

At zero electric field, 6=0, and a weak magnetic field ?C 
which is applied perpendicular to the layer, the orientational 
structure of FN is uniform and planar with orthogonal orien- 
tation of the director and magnetization (n_Lni). With the in- 
crease of the magnetic field this state becomes unstable and 
at a certain critical value the orientational structure of FN be- 
comes deformed, i.e., magnetic Freedericksz transition takes 
place. As shown in Ref.— , for weak segregation {yc > x^) 
magnetic Freedericksz transition is the second-order transi- 
tion, for which He ( |20] i plays role of critical field. At strong 
segregation (x < the orientational transition occurs as the 
first order phase transition (Fig. |5i — curve 1 , Fig. |3}) — point 
Ml). In this case the field //, of equilibrium transition is deter- 
mined by eqn ( |3T1 ). 

Application of electric field £ 1 1 may change the charac- 
ter of orientational transition. The increase of its strength not 
only decreases the transition field He ( |20] |. but also decreases 
the tricritical value >r^ of segregation parameter (l22t . that al- 
lows to satisfy the condition x > >r^. In this case orientational 
transition becomes the transition of the second order (Fig. 
— right curve 2, Fig. |3j) — point M2). For electric field be- 
longs to the interval E^^ <E < E^ + 77 corresponding to am- 
biguity region of the phase diagram (Fig. |3j)), the increase 
of magnetic field leads to a sequence of re-entrant orienta- 
tional phases: nonuniform phase (H < Hi) — uniform phase 
{Hi<H < H2) — nonuniform phase (H > Hj). For E > E^ 
at H = uniform texture of FN becomes unstable and elec- 
tric Freedericksz transition takes place. Application of a mag- 
netic field decreases the deviation of a system from uniform 
orientational state (decrease of director rotation angle inside 
the layer). Subsequent increase of H returns the orientational 
structure of FN atH — Hi into initial uniform state. This effect 
is due to a competition between quadrupole eoea(n • £)^] 
mechanism of the electric field influence on NLC-matrix and 
the dipole jioMsfm ■ IK] influence of magnetic field on 
magnetic moments of ferroparticles. When magnetic field 
is applied, the particles tend to align along the direction of 
the field, but since they are coupled f /(n ' m)^] to NLC- 
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Fig. 5 Director orientation angle (pQ in the middle of FN layer as a 
function of magnetic field H for different values of (a) electric field 
E (Hi = 1.25, H2 = 1.90, H, = 4.72) and (b) segregation parameter 
>i (Hi = 0.67, H2 = 2.59, H, = 1.53); dashed lines represent 
metastable and unstable states 



matrix, they return the NLC-matrix to the initial (uniform) 
state. This state, however, becomes unstable again at mag- 
netic field H = H2, and Freedericksz transition occurs once 
again. In this case magnetic quadrupole jioxd^ ' ^)^] 
mechanism influencing the liquid-crystal matrix prevails over 
the magnetic dipole jioMsfm ■ IK] mechanism. Further in- 
crease of magnetic field gives rise to the increase of FN direc- 
tor rotation angle. 

FigurelSj" represents the director rotation angle ^ in a cen- 
ter of FN layer as a function of magnetic field H for differ- 
ent values of the segregation parameter k at fixed E. Solid 
lines represent stable solutions corresponding to minimal free 
energy of FN, dashed lines — metastable and unstable solu- 
tions. For K > (Fig. I5J) — curve 1, Fig. [3J) — point Pi) 
as the magnetic field H increases the transition from nonuni- 
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form state to uniform at H = Hi and subsequent transition to 
nonuniform state at H ~ H2 occur continuously, i.e., they are 
the second-order phase transitions. With a decrease of seg- 
regation parameter {x < Kjj, Fig. — curve 2, Fig. — 
point P2) the transition from nonuniform state to uniform at 
H = Hi remains the second order, while subsequent orienta- 
tional transition to nonuniform state has a jump at H — Ht dSTT i. 
indicative of the first-order transition. Under further decrease 
of segregation parameter (>f < xr^, Fig.|5j5 — curve 3, Fig.[3j5 
— point P3) at a given value of electric field the re-entrant 
transitions are absent. 

4.2 Electric Freedericksz transition in a magnetic field 

At zero magnetic field, H = 0, the embedded magnetic parti- 
cles behave themselves as a passive admixture. In this case the 
electric field E applied to FN affects only the NLC-matrix, that 
is why electric orientational transition in FN (Fig.|6] curve 1) 
coincides with classical second-order Freedericksz transition 
in 
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Fig. 6 Director orientation angle (po in the middle of FN layer as a 
function of electric field E; Stable branches of solutions are 
represented by solid lines, whereas metastable and unstable — by 
dashed lines 



Presence of a magnetic field || £ can lead to a change 
of the critical field of electric Freedericksz transition and to 
the change of transition character from the second order to the 
first. In weak magnetic fields for x> (Fig.|4|i the transition 
field increases but transition character remains unchanged. 
The growth of a magnetic field up to // = //,„ (where //„, is 
the magnetic field corresponding to the maximum of electric 
field E,„ — E^ + T] on phase diagram Fig. |2]i increases the 
transition field for x < (see Fig. [Hi this transition becomes 
the transition of the first order (Fig. |6] curve 3). With further 
increase of the magnetic field the character of electric Freeder- 
icksz transition in FN does not change but the value of the crit- 



ical field decreases (Fig.|6] curves 4-6). Such nonmonotonic 
behavior of electric Freedericksz transition in FN under mag- 
netic field accompanied by the change of orientational tran- 
sition character is due to the fact that external magnetic field 
IK II £ affects not only the director n [quadrupole magnetic 
mechanism of influence ~ /io;ta(n • IK)^], but also the mag- 
netization m [dipole mechanism ^ ^M^fva ■ 3-C]. Coupling 
between the director and magnetization [term in eqn ^] 
and their mutual orthogonal orientation lead to the competi- 
tion of these mechanisms. In weak magnetic fields the dipole 
mechanism, responsible for orientation of magnetization, pre- 
vails over the quadrupole one, thus increasing the transition 
field. In strong magnetic fields prevails quadrupole mecha- 
nism, responsible for the director behavior, which decreases 
the critical field. The change of transition character is con- 
nected with the dependence of tricritical value of segregation 
parameter x"^ (|27| i on magnetic field, that allows to make x'^ 
either more, or less xe (Fig.|4|i. 

4.3 Ferronematic magnetization 

Magnetization of ferronematic M = M^fva is determined by 
the combination of concentrational / and orientational m dis- 
tributions of magnetic particles in a layer Let us define the 
reduced magnetization of ferronematic ^ - 

M = M/(Mj) =gva = [- g{z) sin v/(z), 0, g{z) cos v/(z)], 

which taking into account the concentrational distribution 
takes the form 

M = mgexp I — cosiJ/- — sin^fffl - V/) I. (32) 
y X X J 

Figure [Tjrepresents the components of reduced magnetiza- 
tion M in FN layer for two values of electric field E. Cho- 
sen value of magnetic field H — 5 is greater than the field 
He = 4.91 of Freedericksz transition in ferronematic at zero 
electric field.— In this case, even in the absence of electric 
field, there are orientational distortions and nonuniformity of 
magnetic particles distribution in the FN layer (segregation ef- 
fect). At £■ = and H (Fig.|2hi curve 1) due to increasing 
of concentration of magnetic particles near the boundaries of 
the layer, is maximal at the boundaries of the layer, while 
the depleted of magnetic particles central part of the FN layer 
is practically not magnetized. Due to rigid planar conditions 
of the director coupling to the boundaries of the layer, \Mx\ 
attains its maximum not on the boundaries, but in-between the 
center of the layer and its boundaries (Fig. |7^, curve 1). 

Application of electric field (Fig.|7] curves 2) contributes to 
further distortion of FN orientational structure and decre asing 
of magnetic particles concentration in the center of the layer, 
which leads to further increase of on the layer boundaries 
and to the growth of |M;r| near them. 
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Fig. 7 Reduced magnetization M of ferronematic layer for x = 2, 
// = 5, <T= 10, yt= 1.56, 6= 10 



Formula (bZb enables us to obtain the average magnetiza- 
tion across the layer: 



Jo 



(My) = 0, (M,) = 2 / gcos xj/ ^5l{(p,ij/)d(p, (33) 
Jo 

which is shown in Fig. |8]as a function of H. In the absence 
of magnetic and electric fields the FN orientational texture 
and distribution of magnetic particles in a layer are uniform: 
n = (1,0,0), m = (0,0,1) and g{z) — 1. In this case aver- 
age magnetization (M) has only z-component (M,) = 1. Ap- 
plication of an electric field E > leads to the distortion 
of FN orientational structure (electric Freedericksz transition) 
and, as a consequence, to the change of average values of the 
suspension magnetization components. For numerical simula- 
tions, electric field has been chosen from the re-entrant phase 
transition interval E^ <E< E^ + rj (see Fig.|2]l. As is seen 
from Fig. [8] application of a magnetic field and the increase of 



its strength returns FN into initial orientational state atH = Hi 
with (Mv) — and (M^) — 1. Further increase of the field 
H leads to the second Freedericksz transition from uniform 
to nonuniform phase at H — H2. At first, there is a decrease 
of (Mj.) and {M^) because of the coupling between the direc- 
tor and magnetization, and then, there appears the asymptotic 
convergence of average magnetization values to their values in 
the undisturbed ferronematic: (Mj) ~ and (Mj.) = 1. 
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Fig. 8 Average value of a reduced magnetization (M) of a 
ferronematic as a function of magnetic field H for >«r = 2, £ = 3.8, 
a = 10, yt= 1.56, fc = 10 

Near the transition from uniform to nonuniform phase the 
reduced magnetization (|32] | can be represented as a powers se- 
ries in (p{z) and (//(z) (fT9l ). where (po^l is determined by re- 
lation (|23] |. In the case of re-entrant phase transitions at £ 7^ 
(see Fig.lU, the Freedericksz transition (direct and re-entrant) 
takes place at H equal to Hi and H2. Then the first nonvanish- 
ing terms of expansion near He for averaged over the layer the 
components of reduced magnetization (l33T l have the form 



-SH(P0, 



(34) 
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from which the nonzero components of the FN dipole part of 
magnetic susceptibility tensor xf^ — d {Mi) / dHj can be written 
as follows 



yd = / yn 



(35) 



where Jh is determined by relation (|23) . 
4.4 Optical phase lag 

Distortion of ferronematic orientational structure can be ob- 
served experimentally by measuring optical phase difference 
of ordinary and extraordinary rays after they pass through 
the ferronematic layer. Let us find out how phase difference 
changes with increasing magnetic field using the following 
formula:— 



o = j [neff[Z) -«oJ dz, 



(36) 



where X is the wavelength of the transmitted beam of 
monochromatic light incident normally to the layer surface, 
n^ff{z) is the effective refractive index, defined by the relation 

1 sin^^(z) cos^ (p{z) 



where n„ and rig are the refractive indices of the ordinary and 
the extraordinary rays, respectively. 

Using condition ( fT4] i. reduced phase lag (|36] l in ferrone- 
matic distorted by applied magnetic and electric fields can be 
written as 



5_ 

8o 



+ cos^(p 
1 - ^ cos^ (p + cos^ ^ 



(37) 

here we introduce the parameter ^ — (n^ — n^) /n^ and the 
phase lag in the absence of magnetic and electric fields ~ 
InLirie — «o) /A. In formula (|37] | the angles Vb and dis- 
tribution function g are given by eqns (|9]l, ( fTTj i, and ( fTSl l, 
whereas the function 'Ji{(p,\j/) — by relation (fT6] l. The re- 
duced phase lag 5/5o (|37| | in ferronematic layer as a function 
of magnetic field H is shown in Fig. |9] 

Note, that phase lag at £ = has been studied previously 
m Ref.5. In the presence of an electric field for weak mag- 
netic segregation (k > k^. Fig. |9] curve 1) a sequence of 
re-entrant magnetic transitions takes place: nonuniform phase 
(H < Hi) — uniform phase (Hi < H < H2) — nonuniform 
phase {H > Hi). The phase lag varies continuously with the 
increase of magnetic field — both transitions are of the sec- 
ond order In the interval H\ < H < H2 ferronematic layer 
has maximal birefringence. Phase lag tends asymptotically 
to zero al H ^ Hjf because of rigid couphng of the director 




magnetic field, H/H^ 



Fig. 9 Reduced phase lag 8 /8o 'm ferronematic layer as a function 
of magnetic field for £ = 3.8, a = 10, /t= 1.56, fc= 10, t, =0.2. 
Here = 1.30, Hi = 0.67, H2 = 2.59, H, = 1.54; solid lines 



correspond to stable states, whereas dashed lines - 
unstable 



' metastable and 



with the layer boundaries. With the decrease of segregation 
parameter {>t < x^, Fig.|9] curve 2) phase lag decreases dis- 
continuously at H = Hi (the first-order transition), and then 
decreases with increasing field. The transition from nonuni- 
form to uniform state at H — Hi is still the transition of the 
second order For smaller values of segregation parameter 
(x ^ x^i, Fig. |9l curve 3) there is no re-entrant orientational 
transitions. Application of a magnetic field initially increases 
the phase lag, which never reaches maximum value d ~ 5q, 
because FN orientational structure does not return to initial 
uniform state. Further increase of magnetic field leads to the 
decrease of phase lag. 

5 Conclusion 

In the present paper we have studied the influence of mag- 
netic and electric fields on a suspension of magnetic parti- 
cles in nematic liquid crystal (ferronematic). Both fields are 
able to induce the orientational Freedericksz transitions in fer- 
ronematics independently of one another However, the re- 
sponse of a suspension on these fields is different. An ex- 
ternal magnetic field affects both the orientation of magnetic 
particles (dipole mechanism) and the liquid crystal structure 
(quadrupole mechanism). Coupling the director with mag- 
netization leads to competition between these mechanisms of 
influence. In contrast to magnetic field an electric field can 
directly affect one of the subsystems — liquid crystal. We 
have shown that in certain electric field range the magnetic 
field can induce a sequence of re-entrant orientational transi- 
tions in ferronematic: nonuniform phase — uniform phase — 
nonuniform phase. This phenomenon is caused by the inter- 
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play between the dipole (ferromagnetic) and quadrupole (di- 
electric and diamagnetic) mechanisms of the field influence on 
a ferronematic structure. We have found that these re-entrant 
Freedericksz transitions exhibit tricritical behavior, i.e., they 
can be of the first or the second order The character of 
the transitions depends on a segregation degree and external 
fields. Moreover, electric field can change the character of the 
magnetic Freedericksz transition in ferronematic from the first 
order to the second, but not vice versa. 

We have shown that re-entrant Freedericksz transitions in 
ferronematic can not be induced by electric field at any fixed 
values of a magnetic field, since the threshold electric field 
for the orientational transition is the single-valued function 
of a magnetic field. The presence of a magnetic field can 
change the type of electric Freedericksz transition in ferrone- 
matic from the second order to the first only. In addition, we 
reveal a nonmonotonic behavior of the threshold electric field 
for this equilibrium transition. 
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